Simulating spin-charge separation with light 
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In this work we show that stationary light-matter excitations generated inside a hollow one- 
dimensional waveguide filled with atoms, can be made to generate a photonic two-component Lieb 
Liniger model. We explain how to prepare and drive the atomic system to a strongly interacting 
regime where spin-charge separation could be possible. We then proceed by explaining how to 
measure the corresponding effective spin and charge densities and velocities through standard optical 
methods based in measuring dynamically the emitted photon intensities or by analyzing the photon 
spectrum. The relevant interactions exhibit the necessary tunability both to generate and efficiently 
observe spin charge separation with current technology. 

PACS numbers: 



One of the most counterintuitive characteristics of one 
dimensional electron gases is spin-charge separation. In 
this case the electrons cease to behave as single particles 
comprised of spin and charge. Instead collective excita- 
tions appear carrying only charge (and no spin) or only 
spin (and no charge) which can propagate through the 
system with different velocities [1. Spin charge separa- 
tion was predicted in tunneling experiments in metal- 
lic chains [2 j, organic conductors [3], carbon nanotubes 
[4] and more recently in quantum wires [5 . However 
due to the complexity of the structures used, measur- 
ing the spectral function and observing distinct spinon 
and holon branches has yet to be conclusive. In parallel 
with these works, artificially engineered many-body sys- 
tems that could simulate condensed matter effects in well 
controllable environments have been developed the last 
decade. Cold atoms and ions traps are the most famous 
example [6-8 , and strongly interacting photons (SIP) 
the most recent developments. Proposals to observe spin 
charge separation have also been in place for some time in 
cold atoms, including both bosonic and fermionic species 
[T5H20] . However the lack of necessary individual acces- 
sibility and measurement, and the challenges in trapping 
and cooling especially fermionic gases make current re- 
sults inconclusive so far 

SIPs on the other hand, as hybrid light-matter quan- 
tum simulators promise to provide the necessary extra 
manipulation and measurement lacking so far from other 
proposals. The efficient quantum simulation of photonic 
and polaritonic Mott transitions and the crystallization 
of photons was shown to be possible using both arrays 
of coupled resonators or stationary polaritons in atomic 
gases [9HT4]. We show here for the first time that spin- 
charge separation could be efficiently observed in a such 
strongly correlated quantum optical system. 



In parallel with the seminal works in Luttinger liq- 
uid and in cold atoms physics, significant progress has 
also been made in a different field. Quantum nonlinear 
optics has shown that strong interactions between light 
pulses, even at the single photon level is possible with 
numerous applications in photon switching and quan- 
tum memories. Employing Electromagnetically Induced 
Transparency techniques (EIT), few photons pulses can 
be slowed down, stored in atomic medium and then ef- 
ficiently recovered again in an all coherent manner [2T1 - 
[23] . In this work we show that stationary light-matter 
excitations generated inside a hollow one-dimensional 
waveguide filled with atoms, can be made to generate a 
photonic two-component Lieb Liniger model. Moreover 
the relevant interactions exhibit the necessary tunabil- 
ity both to generate and efficiently observe spin charge 
separation. 

Consider a waveguide, a hollow fiber for example filled 
with two types of atomic gases a and b. Two quan- 
tum light fields and two classical fields Ei t ±,E2,± f2i,± 
and 0,2, ± respectively can propagate towards the left and 
right directions and couple to the atoms as shown in 
Fig. 1. Initially two resonant optical pulses carried by 
^i )+ ,£^2,+ are incident from one direction, say the left 
side. They are injected into the waveguide with the co- 
propagating control fields and ^2,+ initially turned 
on. The Hamiltonian in the Schrodinger picture can be 
expressed as a sum of two independent parts describing 
the evolution of the different atomic species H = H a + H b 
with 
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with x = {a, b} labeling the two atomic species a and 
b. Here, the continuous collective atomic spin operators 
a pq ~ ^pq ( z it) (PiQ = {I? describe the average of 

|p) x (q\* over the atoms of type x in a small but macro- 
scopic region around spacial coordinate z. The densi- 
ties of the two species atoms in the same region are as- 
sumed to be different and equal to n a and n b respectively. 
Moreover E^±(z,t) and £\±(t), with i = {1,2}, are the 
two quantum and two classical fields slowly varying op- 
erators with frequencies cjqu and and wave vectors 



kqu and k c ^ 
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respectively. Both quantum fields E\ t ± and 
£^2,± drive four possible atomic transitions and for sim- 
plicity we assume that each field couples with the same 
strength to the two transitions of the different atoms. 
The field Ei y ±(z, t) is detuned by A| from the transition 
|l) a |2) a and by A| from |3) a |4) a . Similarly, the 
quantum field E2,±{z,t) is off-resonant from 



^|2) b 



and |3) b -> |4) b by A b and A b as shown on F 
Finally, the applied classical control beams with Rabi 
frequency fii } ±(t) and ^2,±(^) couple to both atoms and 
drive the transitions |3) x — >• |2) x . The evolution of the 
slowly- varying quantum operators E iy ±(z,t) is given by 
four Maxwell-Bloch equations 




dt dz 



and 



Ei,±(z,t) = \\Fhtn%g\ (^ 2}± (z, t) + <rf 4 

(2) 
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where we have introduced the slowly varying collec- 
tive operators a^ q 



34, ± 



and at = a b ^U, *)e ife S )z + at _ fle" 1 ^*. In the 



derivation of Eqs. ( |2|3[ ) we have ignored the high- 
frequency terms oscillating at frequency 2A x t. The 
Langevin-Bloch equations which define the evolution of 
the spin operators a* q can be derived following the 
methods in [23 We proceed to define the polariton 
operators as S&i^± = cos6 a Ei^± — sin a y /27rn^a^ 1 and 
cos# b I?2,± — sin h ^2i\n° z a\ x where tan# a = 

For sim- 



FIG. 1: A schematic of system under consideration. A hol- 
low fiber with two quantum light fields , ^2,+ , and two 
pairs of classical fields Qi,± and Q2,± propagating towards 
the left and right directions. The fiber is filled with two 
atomic gases a and b. Initially two coherent few photon op- 
/ tkal pulses carried by Ei t +,E2,+ are incident from one di- 
^ recxi'on, say the left side. By adiabatically switching off the 
corresponding control fields the quantum pulses can be slowed 
down and eventually mapped to stationary light-matter ex- 
citation (polaritons). Appropriate tuning of the couplings of 
the light fields to the corresponding atomic transitions, forces 
the trapped polaritons to behave as ID quantum liquid of 
{z\^pj types of bosons obeying the Lieb Liniger modelpQ. The 
available tunability of the effective interaction parameters al- 
lows us to reach the spin charge separation regime. Subse- 
quent release of the trapped polaritons through their coherent 
mapping back to propagating light pulses, allows for the effi- 
cient measurement of both the dynamics of the propagation of 
the spin and charge quasiparticles or of the spectral function 
characteristic of the effect taking place. This is done though 
standard quantum optical measurement techniques measur- 
ing correlations in the intensities of the output fields and by 
analyzing the optical spectrum 
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plicity we have assumed that the amplitudes of the 
counter propagating classical fields are equal, i.e., f^ ? ± = 
= 1,2). In the limit when the excitations are 
mostly in the spin- wave form, i.e. sin# x ~ 1, and since 
^31 = -5iA,±M and cr^ = -g%E 2 ,±/p2, the po- 
lariton operators equal &i,± = g\^hrn^E\^±jVt\ and 

Setting ^i )2 = Oi, 2 ;++^i,2;-)/2 andAi )2 = Oi,2 ; +- 
^ r i j 2 ; -)/2 as the symmetric and antisymmetric combina- 
tions of the two polaritons, we derive the equations of 



motion for the polariton combination ^i, A\ : 
d t *i + vV>d x A l = - v^^n^tfi 
, (5f) 2 
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A similar equation holds of the pair ^-2- Vg 1 ^ = 
*/«Q2/7r(<tf) 2 < and ^ 2) = v^^\l^g\fn h z are the 
corresponding group velocities of the propagating po- 
laritons and z/ 1 ' 2 ) are the velocities for each quantum 
field in an empty waveguide. The noise terms in Eqs. 
Q account for the dissipative processes that take place 
during the evolution. However for the dark state po- 
laritons under consideration for the case of degener- 
ate lower levels |1) and |3), and T « \A%\, \A%\ the 
losses in the timescales of interest are not significant 
and thus can be neglected [12] [24] . Assuming sufficient 
optical depth of a few thousand and a large ratio be- 
tween the number of atoms to the number of photon 
nl/po.i ~ 10 4 , the antisymmetric combinations A\ and 
A2 can be adiabatically eliminated from the equations 
of motion for the polaritons and moreover, the nonlinear 
terms ^\^f 1 A 1 , ^\^! 2 M are negligible. Then, Eqs. Q 
simplify to d t V 1 = 2^7 9^ 1 + U 1 ^\^ 1 ^ 1 + Vi^^^i 

and d t ^2 = ^^^2 + U 2 ^l^2^2 + V 2 ^\^i^2 These 
are the equation of motion of a two-component Luttinger 
liquid of polaritons with the corresponding Hamiltonian: 
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The interaction parameters relating to the effective ki- 
netic energy, intra and inter species repulsion can be 
tuned by controlling the external control fields and the 
various detunings. ra$ is effective mass for the ith po- 

lariton with — 



and ± 



- r, h K with 



1D the spontaneous emission rate of a single atom in 
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the waveguide modes. The intra repulsion are given by 

r 1D ^> ^ r 1D „< 2 > 
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and the inter-species by V\ 2 — 

(^) 2 A^(2) 



with V 2 
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Let us n ow define th e effective pa rameters ui = 
yJp^iUi/rrii and K { = V / ^P(W ( m iUi) with p ,i equal 
to n p h,i, the photon number of ith quantum field in our 
scheme. It is known from the works by Girardeau and 
others, that when u\ = u 2 = u and K { = K 2 = K the 
above LL Hamiltonian Eq. 1 can be transformed into new 
one comprised separately of two parts, the charge part H c 
and the spin part H s as H = H c + H s [1 . This allows 
for the separation of a single excitation into two separate 
ones each comprised of spins or of charge/density. These 
can propagate through the liquid with different veloci- 



ties given by u CjS 



u\l± 



(Vi+V2)K 



In our case, the 



charge (spin) density corresponds to the sum(difference) 
of the corresponding polaritons densities n c ^ s = n\ ± ri2 
with m = and ^ = (^ ;+ + *i ; -)/2 the sym- 

metric combinations of the two counter-propagating dark 
state polaritons generated by each set of atom-field in- 
teraction. 

We note here that in our case, Eq. 1 describes a non- 
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FIG. 2: The single particle spectral function for our photonic 
system. The effective spin and charge velocities are u s = 0.5 
and u c = 1 and can be achieved for optical depths around 
OD=4000, with 10 photons in each pulse initially, and single 
atom cooperativity for each atomic species of 0.2. Subsequent 
releasing of one of the trapped polaritons through its coherent 
mapping to a propagating light pulse, allows for the efficient 
measurement of the spectral function. This is done though 
standard quantum optical measurement techniques measur- 
ing correlations in the intensities of the output fields and by 
analyzing the optical spectrum. In the inset we plot a 2D cut 

for the value of q=2. The units of u and q are in x 2^j^u 

and tt/zq respectively. 



equilibrium situation where we follow the dynamic evo- 
lution of two polariton pulses. Here the effective masses 
mi 5 2 are kept constant in time while the effective inter 
and intra species repulsions and V12 are increased to 
reach the necessary regime. This is possible by keeping 
for example the corresponding detunings A§, A2 con- 
stant, while shifting the level 4 responsible for the non- 
linear shifts by changing the detunings A^A^. Labeling 
71 5 2 the corresponding ratio of the interaction to the ki- 
netic energies for each polariton species, at the time of 
the pulses entering the fiber, conditions are such that 
71 5 2 < 1. At this phase interactions are weak and the 
photons expand freely due to dispersion. By shifting in 
the 4th- levels, the system can be driven into the strongly 
interacting regime 71,2 > 1 necessary for spin charge sep- 
aration to occur. [26] . 

We will now analyze in detail the steps characteriz- 
ing the preparation, release, and the measurement of the 
separation, estimating the quantum optical parameters 
required to reach the LL regime. First the fiber is in- 
jected with the two quantum coherent pulses £"2,+ 
from the left and the adiabatic switching off the two clas- 
sical fields and ^2,+ occurs. The polaritons ^1 and 
^2 are formed and get trapped. Then, all four classical 
fields f2i } ± and ^2,± are slowly switched on so that an 
effective Bragg grating is formed. The latter forces the 
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polaritons - which now have acquired a small photonic 
component as well to evolve under the LL model. The 
interaction parameters (couplings and detunings) char- 
acterizing the intra- and inter- species repulsion can be 
tuned in time in order to reach the necessary regime for 
density separation. Once this is achieved one field, say 
say f2i+, from the pair of control fields that trap po- 
lariton type is slowly turned off. This will release 
the corresponding quasiparticle and allow it to propa- 
gate towards the exit of the fiber. As all correlations es- 
tablished in the previous step- the evolution under the 
LL Hamiltonian Eq. 1-are retained, this wavepacket 
containing mostly of light now, comprises of two sepa- 
rated parts propagating with different velocities u s and 
u c towards the end of the fiber. There are two ways 
to observe the separation in this system. The first is 
by measuring the time evolution of a single excitation 
as usually predicted for the case of cold atoms propos- 
als as well[15] [19j [20]. In our case the charge(spin) 
density waves, after the release of the polaritons, will 
transfer to the sum(difference) of the corresponding time 
dependent photon intensities for each propagating field 
^c,s = ni ± ri2 with rii = (£^ + £^ ;+ ). Using standard 
optical measurement techniques one could dynamically 
probe the time evolution of the corresponding effective 
charge and spin densities maxima and from that infer 
the corresponding velocities. The second way to measure 
the effect is through directly probing the correlations es- 
tablished during the evolution with the LL Hamiltonian 
by measuring the corresponding single particle spectral 
function as usually proposed in condensed matter exper- 
iments S(q,uj) p~H5]. In our case, as the quasiparticle 
after the coherent release are propagating photo pulses, 
this is directly accessible by measuring the first order cor- 
relation function ( which corresponds to analyzing the 
spectrum) of one of the initial fields, say E\ for a spe- 
cific quasimomenta q. This should exhibit two peaks as 
shown in the inset of Fig. 2 centered around u s q and u c q 
with q the corresponding quasi-momentum of the propa- 
gating excitation. In our case q is inversely proportional 
to the initial extend of the pulse zq. For clear distiction 
between the two spin and charge peaks we should tune 
our optical detectors around q = 2tt/zo . In Fig. 2 we 
plot the spectral function as analytically derived for a 
two component bosonic system [25]. for intra and inter 
species values corresponding in our system to u = 1 and 
k = 1. Here U\/V\2 — U2/V12 = 0.6 which in turn give a 



ratio between the effective charge and spin velocities of 
u c /u s = 2. The corresponding necessary values of the pa- 
rameters in our system are optical depths OD=4000 and 
single atom cooperativity of 0.2 and roughly 10 photons 
initially each pulse [27]. We note here in calculating the 
latter we took into account that both the linear and the 
nonlinear loss mechanisms will define a maximum evolu- 
tion time t max for the second phase of the process. This 
in turn will give the following condition on the achievable 
ratio of interaction to kinetic energies [12] [24] 

12 . / ,P\^2 h \ Q r OD a >\ , . 

7max ~min(7 exp(— —,77/?— y—f^)i (6) 

1^2 I 7V ph 

Optimizing over A2' 6 will give the numbers mentioned 
above. We add here that in order to omit the noise terms 
and also be able to neglect the higher order derivatives 
so that the system is within good approximation a two 
component Luttinger liquid as described in Eq. 1, the 
following conditions need to be satisfied. The Stark shift 
of levels 1 2) for each atomic species should stay within 
the EIT transparency window and a bound on spin-wave 
excitations wavevector should be in place. These con- 
ditions translate to n\/n% < |A||/|r - 2iAg|,n 2 /n* <C 
|Ai|/|r - 2iA|| and n x jn a z <C T 1D /\2A% + iTln 2 /n b z <C 
T/|2A2 + iT\ which is satisfied in our case as the num- 
bers of atoms can be two or three order of magnitude 
larger than the number of photons. We note that in this 
case the ratio between kinetic and repulsion energies for 
each species can reach the value of 40 which for the sin- 
gle component case was shown to lead to a Tonks gas 
of photons [12]. For spin-charge separation, the required 
repulsions could be of a smaller value, thus relaxing the 
overall quantum optical requirements but the two peaks 
in the spectral function will be less pronounced. 

In conclusion we have shown that stationary light- 
matter excitations generated inside a hollow one- 
dimensional waveguide filled with atoms, can be made to 
simulate a photonic two-component Lieb Liniger model. 
Moreover the relevant interactions exhibit the necessary 
tunability both to generate and efficiently observe spin 
charge separation using standard quantum optical meth- 
ods. 
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